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Solid substrates can be endued with self-organized regular stripe patterns of nanoscopic lengthscale 
by Langmuir-Blodgett transfer of organic monolayers. Here we consider the effect of periodically 
prepatterned substrates on this process of pattern formation. It leads to a time periodic forcing of 
the oscillatory behavior at the meniscus. Utilizing higher order synchronization with this forcing, 
complex periodic patterns of predefined wavelength can be created. The dependence of the synchro- 
nization on the amplitude and the wavelength of the wetting contrast is investigated in one and two 
spatial dimensions and the resulting patterns are discussed. Furthermore, the effect of prepatterned 
substrates on the pattern selection process is investigated. 



I. INTRODUCTION 

Prepatterned substrates are widely used as means to 
control wetting and dewetting processes PHI]- Chemi- 
cally or topographically prestructured, the substrate ex- 
hibits heterogeneous wetting properties and thereby en- 
forces an extrinsic pattern upon a covering liquid layer. 
A particularly interesting scenario emerges if the liq- 
uid itself exhibits pattern formation intrinsically, yield- 
ing regular or irregular structures even on homogeneous 
substrates. This is the case, for example, in Langmuir- 
Blodgett transfer of organic and metallic monolayers [5j 
[6] and in controlled evaporation of nanoparticle suspen- 
sions [7 j. On a prepatterned substrate, this gives rise to 
the formation of a final pattern resulting from a com- 
petition between the extrinsic structure imposed by the 
substrate and the intrinsic structure that is formed nat- 
urally in the system. 

Here, we want to focus on the effect of periodically 
prepatterned substrates on the transfer of monolayers of 
amphiphilic molecules like pulmonary surfactant DPPC 
(dipalmitoylphosphatidylcholine) . Using homogeneous 
substrates, this transfer process results in the sponta- 
neous formation of nanochannel lattices, that is, regu- 
lar stripe patterns with a typical wavelength of merely 
a few hundred nanometers [SJ [HJ |9]. The correspond- 
ing experiments utilize the Langmuir-Blodgett method, 
where the substrate is withdrawn from a trough filled 
with water on which a monolayer has been prepared, 
and are conducted under conditions close to the first or- 
der phase transition between the liquid-expanded (LE) 
and the liquid-condensed (LC) phases of the monolayer. 
The observed patterns consist of ordered arrays of LE 
and LC domains, arranged into stripes oriented either 
parallel or perpendicular to the contact line. For certain 
experimental parameters, one can even obtain a combina- 
tion of both orientations, resulting in a rectangular struc- 
ture. The formation of the nanochannels results from 
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a substrate-monolayer interaction, commonly referred to 
as substrate- mediated condensation (SMC) [9HT3]. and 
can be understood theoretically as a phase decomposition 
process under influence of the corresponding interaction 
field M]. 

Since the menisucs, where the condensated domains 
are generated, moves across the prepatterned substrate, 
the spatially periodic prestructure is translated into a 
spatiotemporal forcing. The effect of spatiotemporal 
forcing on different pattern forming systems, e.g., photo- 
sensitive Turing systems [15], has been investigated ear- 
lier. However, in these cases, an external frequency is 
imposed in an otherwise nonoscillatory system, whereas 
even the unforced pattern formation in monolayer trans- 
fer is an oscillatory process. 

The theoretical investigation of dewetting systems with 
surfactants is typically based on the lubrication approx- 
imation which greatly facilitates the description of thin 
film flow [T6HT8] . In this framework, surfactant-covered 
thin liquid films can be modeled by two coupled non- 
linear differential equations describing the height profile 
of the film and the surfactant density [T7H20] . whereas 
prepatterned substrates can be described by a spatially 
varying disjoining potential [2TJ[22]. In addition, mono- 
layers in the vicinity of the LE-LC phase transition can 
be modeled by coupling of the lateral pressure, that is the 
surface tension, to a Cahn-Hilliard type free-energy func- 
tional [23]. The substrate- mediated condensation can 
then be described by an interaction field acting on the 
monolayer [T4j . whose strength depends on the distance 
to the substrate, that is, the thickness of the liquid film 
dividing substrate and monolayer. Due to this interac- 
tion, the condensed state is energetically favored when 
the monolayer is close to the substrate, leading to the 
partial condensation necessary to form the alternating 
LE/LC stripes. 

Investigating the influence of periodically prepatterned 
substrates on this pattern formation process we shall 
find synchronization phenomena [24 of various orders, 
leading to the formation of complex periodic patterns 
depending on the amplitude and the wavelength of the 
substrate's wetting contrast. We will further show that a 
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FIG. 1. Schematic drawing of a surfactant- covered meniscus 
on a prepatterned substrate. The height profile h(x, y, t) indi- 
cates the film thickness at location (x, y) and time t, whereas 
7(x, y, t) describes the surfactant density at the surface above 
(x,y). 



periodic prepattern can be used to inhibit an orientation 
transition from stripes arranged parallel to the contact 
line to perpendicular stripes and even allows to produce 
new structures such as well-ordered domains of circular 
LE domains alternating with LE/LC stripes. 



II. GOVERNING EQUATIONS 

A. Thin film flow 

Let us consider an evaporating thin film of water on a 
periodically prestructured substrate which is withdrawn 
with constant velocity v in negative ^-direction from a 
water reservoir which ensures a constant film height at 
the right boundary of the system (see figure [I]). A menis- 
cus is then given by the balance of evaporation and the 
supply of fresh water, which is carried from the reservoir 
by the moving plate, yielding a contact line moving rel- 
atively to the substrate. The water reservoir is assumed 
to be covered with a surfactant monolayer of constant 
density, which is then carried by the flow towards the 
contact line. 

The liquid film is described by its height profile h(x,t), 
which indicates the local film thickness, whereas the sur- 
factant density at the surface above the point x — (x, y) 
is described by the function j(x, t). Introducing the 
scales ho, Zo, and 70 for height, length, time, and 
surfactant density, the time evolution equations for the 
height profile H(X,T) := h(x/lo,t/to)/ho and the sur- 
factant density T(X,T) := j(x/Iq : £/£o)/7o can be writ- 
ten in dimensionless form as [14] 
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with the generalized pressure 

and E = eCa -1 <5\ Here, the parameter e = Hq/Iq signifies 
the ratio of the characteristic height- and length scales of 



the system and is assumed to be small in thin film prob- 
lems. The inverse capillary number Ca _1 = cro/(rjUo) is 
the ratio of a characteristic surface tension ctq, the dy- 
namic viscosity 77 and the characteristic velocity Uo = 
lo/to. The scaled surface tension a = cr/ao of the water 
film depends on the density T of the surfactant mono- 
layer. Therefore, spatial variations of T lead to spatial 
variations of <r, inducing Marangoni forces which are ac- 
counted for by the terms ~ VS in the equations and 
([2]). The dependence of a on the monolayer density can 
be described in terms of a lateral pressure [26] and, as 
will be outlined in the next section, is deeply related to 
the surfactant thermodynamics. 

Since the substrate is withdrawn with a constant ve- 
locity both equations and Q include an advective 
term involving V = v/Uq. 

Note, that the generalized pressure P comprises not 
only the Laplace pressure ~ V 2 i7, which is present in 
every free fluid interface problem, but also the disjoining 
pressure II (H) due to the interaction of substrate and 
liquid. These interactions become significant for liquid 
films of heights lower than a few hundred nanometers. In 
the literature, different expressions for H(H) have been 
considered (see [18] and references therein for a discussion 
of possible choices). Here, we employ the dimensionless 
expression 
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Hamaker constants as^a^ > 0, yielding, for M = const., 
a precursor height defined by H(H p ) = <^> H p = 
yj Aq/As. The function M(X) accounts for the spatially 
heterogeneous wettability of the prepatterned substrate. 
This approach is similar to the one in [2TJ[22]. Since we 
are interested in periodically prestructured substrates, we 
choose M(X) to be a periodic modulation 

M(X) = l + p£(X + FT) 

which moves along with the substrate and has the form 
of a kink-antikink train 
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with the periodic triangle function 
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where {x} := x mod 1 denotes the fractional part of 
x. This means, that we have a modulation with ampli- 
tude p and wavelength L p . For p — >• the homogeneous 
substrate is obtained. The constant ft p determines the 
sharpness of the contrast. It has to be noted, that very 
small ftp would violate the assumption of slow variation 
in x- and ^-direction, which is the basis of the lubrication 
approximation. 
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Due to evaporation, a sink term is present in the evo- 
lution equation for H, so that the fluid volume is not 
conserved. Here 5ji = /i w — /i v denotes the difference of 
the chemical potentials of the water film and the ambient 
vapor phase, whereas Ev = rjl^Qe/hQ is the evaporation 
number with effective rate constant Q e . The pressure in 
the vapor above the film is assumed to be close to the 
saturation pressure, allowing us to identify the chemical 
potential of the water film with the negative generalized 
pressure [27l[28], that is /i w = — P, fi Y = const. 

It has to be noted that monolayer rigidity is not taken 
into account in the model given by eqs. ([!]) and (2|. This 
is due to the fact, that, going from low to high densities, 
a significant surface rigidity K c is observed experimen- 
tally only after the LC to solid (S) phase transition. To 
give a more precise estimation of the influence of surface 
rigidity, one can supplement lubrication theory by a con- 
tribution resulting from minimization of the correspond- 
ing Helfrich energy of a stiff interface. As a result, an ad- 
dtional quantity a 2 = K c /Iq enters the equations with ex- 
actly the same prefactors as the surface tension a. Thus, 
(J2 can be regarded as a second surface tension and can 
be directly compared in its magnitude to a. Even in the 
S phase, which is significantly more rigid than the liquid- 
condensed phase considered in our problem, the bending 
rigidity of a DPPC monolayer is only K c < 70 k^T [25] , 
that is of order 10 -19 J. Since lo ~ 10 _T m one can es- 
timate that in this case 02 ~ 10 -5 J/m 2 , which is three 
orders of magnitude below the surface tension of water 
and therefore neglect able. For other materials, however, 
the rigidity can be significantly larger, so that it might 
have to be taken into account for an accurate description 
of thin films covered with such monolayers. That is, the 
applicability of the description presented here is limited 
to materials where the rigidity K c of the monolayer is 
less than a few hundred k^T. 



for a general non- homogeneous system [29 , i.e. ^[r] = 
J d 2 XF(T(X)) with a free-energy density 
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Here, K = / {KqCFq) denotes the nondimensionalized 
form of the line tension n between domains of different 
densities, while the function Fh om represents the free- 
energy of a system with homogeneous surfactant density 
T. To focus on the LE-LC first-order phase transition 
of the monolayer, we consider a bistable system where 
the low- and high-density minima of Fh om correspond 
to the liquid-expanded and the liquid-condensed phase, 
respectively. For the sake of simplicity we assume the 
function Fh om to be a symmetric double-well, centered 
around a critical density r cr . It has to be emphasized, 
that this assumption is no real limitation, since, around 
the coexistence plateau, a wide range of experimentally 
obtained pressure-area isotherms can be fitted under the 
assumption of a symmetric Fh om - Using T = T — r cr , we 
write 
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Here, Mi equals the curvature of the function Fh om in 
its center T = and M2 denotes the difference between 
the densities FlejIYc of the pure LE and LC phases and 
the critical density r cr , that is M2 := |Tle — r cr | = 
ITlc — r cr |. The constant F cr is the free-energy of a 
homogeneous monolayer with the critical density r cr . 

The lateral pressure can be obtained from the surfac- 
tant equation of state, which is related to the free-energy 

by 
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B. Surfactant thermodynamics and substrate 
mediated condensation 

The flow of a thin liquid film on a substrate depends on 
the surface tension a at the free liquid- vapor interface. 
In the presence of a surfactant monolayer, the surface 
tension varies across the surface, depending on the local 
surfactant density T and its spatial derivatives. This de- 
pendence is described by the lateral pressure P\ at and the 
relation [26l 
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where a a b s = crabs/00 denotes the scaled surface tension 
in absence of any surfactant. 

The relation between lateral pressure and surfactant 
density is determined by an equation of state which can 
be obtained from the free-energy of the monolayer. Since 
the density V may vary across the surface, its free-energy 
can be written in the form proposed by Cahn and Hilliard 



where the chemical potential /i is obtained from T by 
functional derivation with respect to T, yielding 
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Inserting eqs. ^ and ^ into ([5|, the physical meaning 
of F cr becomes obvious: 

^cr = c(r c r) — ^abs- 

The above description is still incomplete as inhomoge- 
neous systems are considered. In case of a spatially vary- 
ing density T the pressure is given by a tensor P with 
the components [30j [31] 



Pij =Pi a tS ij +e 2 K(d i T)(d j T). 



(10) 



However, in the Laplace pressure, which is of order e 3 , 
the derivative terms of P are of fifth order in e. Since 
e is assumed to be small, they can be safely ignored. In 
the Marangoni term VE, the lateral pressure enters in 
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the form V • P = V (P ]at + e 2 if(Vr) 2 
can simply use the scalar pressure 
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Therefore, one 
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throughout all calculations. 

It is an experimentally well-established fact that inter- 
action with the substrate influences the thermodynamic 
behavior of the monolayer, when the liquid subphase di- 
viding them becomes very thin [9HT3]. It is by this in- 
fluence, that a pure monolayer condensates at the sub- 
strate even under conditions where it would still be per- 
fectly stable on a subphase of ~ 100 nm thickness. This 
substrate-mediated condensation effect can be modeled 
by inclusion of an external field in the free-energy of the 
monolayer, whose strength depends on the distance to 
the substrate which, in our case, is given by the thickness 
of the water film H. The free-energy density in presence 
of the substrate-monolayer interaction is then given by 
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Here, S(H) describes the functional dependence of the 
strength of the substrate-mediated condensation on the 
distance. Unfortunately, the exact form of S(H) still 
remains to be determined by experiment or from a mi- 
croscopic theory. Nevertheless, there exist rough esti- 
mates of the overall lateral pressure difference caused by 
SMC: For DPPC and DMPE the LE-LC coexistence pres- 
sure at the substrate can be several millinewton per me- 
ter lower than for a floating monolayer [T2] . yielding an 
overall difference of more than 40%. In terms of S(H) 
this means that S(H) goes to a finite negative value for 
H ^ H p , leading to a tilt of F\ lom towards its higher 
density minimum, that is the LC phase, as the mono- 
layer approaches the substrate. Also, S(H) should vanish 
quickly for H — >• oo, since the SMC can only be observed 
very close to the substrate. 

Although the results presented in this article are quali- 
tatively valid for a class of functions S(H) meeting these 
minimal constraints, without loss of generality we choose 
S(H) = B^(H), where ^(H) = J dH'IV(H') with inte- 
gration constant zero is the potential of the substrate- 
liquid interaction and B is a positive coupling constant. 
By this choice, SMC acts on length scales comparable to 
the substrate-liquid interaction. 



time scale of the problem is obtained from the relation 
to = 7o^ 2 77/(/iq(Jq). With this choice of scales, e obtains 
a direct physical meaning as 



1 (Tp 

7o^ 



2/3 



(13) 



relates the forces acting at the monolayer covered surface, 
characterized by 00,70, and ft, to the forces acting be- 
tween substrate and liquid, determined by the constants 
a 3 and a§. 

Since we are interested in systems in vicinity of the 
LE-LC phase transition, it is convenient to set the scale 
7o equal to the critical surfactant density, that is 70 = 7 cr 
and to define the characteristic surface tension <jo = 
c"(7cr) as the surface tension of water covered with a 
monolayer of density 7 cr . The resulting scaled expres- 
sion for the surface tension is given by 
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Analogously, the generalized pressure P in the equations 
([T]) and ([2| is obtained as 

P = (1 - e 2 P hom )V 2 tf - 11(H) (17) 

whereas VS can be written in the form 

vs = -VP hom + e- 2 rv 3 r. (is) 

Note, that P does not contain any derivatives of T, due 
to the argument following equation (10). 



III. 



TRANSFER ONTO PERIODICALLY 
PREPATTERNED SUBSTRATES 

A. Methods 



C. Choice of scales 

So far, we did not specify a certain set of scales ho, lo, 
to, and 7o. A convenient choice results from demanding 
e 3 Ca _1 = e 2 K = H p = 1. This means, that we use the 
precursor height as a natural height scale of the system 
by setting ho = h p . Due to this choice we have A := 
A 3 = A 6 = (loKa3)/((T%h*) in equation 

The characteristic length scale lo of the system is then 
given by the ratio of the surface tension <Jo and the do- 
main line tension k as lo = 70 VW^o and the natural 



The system ([TJl-Q is simulated numerically in one 
and two spatial dimensions using second-order finite dif- 
ferences for discretization of space and an embedded 
adaptive Runge-Kutta scheme of order 4 (5) for time- 
stepping [32] . In the ID (2D) case, the integration do- 
main is [0, L] ([0, L] x [0, L\) with L = 1200, discretized on 
a grid of 384 (384 x 384) gridpoints. For a more detailed 
investigation of some of the obtained patterns, a larger 
domain with L = 2000 is simulated using 640 gridpoints. 
In X-direction the boundary conditions are 

T(L) = T L , d 2 x T(L) = = d x T(0) = d 2 x T(0), 
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H(L) = H L , d 2 x H(L) = = d x H(0) = d 2 x H(0), 



with = 0.835, corresponding to transfer from a reser- 
voir which is covered with a homogeneous monolayer 
in the pure LE phase. The choice of boundary condi- 
tions at X = is not obvious, since ideally the material 
would simply leave the integration domain to the left. 
However, such nonreflective boundary conditions are ex- 
tremely complicated to realize for equations like and 
([2}. Therefore, we use the above simple set of condi- 
tions which we tested on three different domanins with 
L = 1200, 2400, and 3600 without finding any influence of 
the boundaries on the pattern formation at the meniscus. 
In the two-dimensional simulations, periodic boundary 
conditions hold in the F-direction. 

We investigate periodically prepatterned substrates 
with wavelengths L p = 200, 300, and 400, each for wet- 
ting contrast amplitudes between p = 0.001 and p = 0.01. 
All three prepattern wavelengths are on the same order 
of magnitude as the natural wavelengths of the system. 
For comparison, we will start with a discussion of the 
case of homogeneous substrates, p = 0. 

For each value of p, one can find a range of trans- 
fer velocities V for which the surfactant monolayer is 
deposited in regular patterns of periodically alternating 
stripes of liquid expanded and liquid condensed mate- 
rial. The wavelength of the patterns, as well as the ra- 
tio of the widths of the LE and LC domains depend on 
V. The patterns consist of travelling domains, which 
are created at the contact line and then advected by 
the moving substrate (see figure [5] for a spacetime plot 
of patterns evolving on a homogeneous substrate as ob- 
tained from one dimensional simulations). Away from 
the contact line, the surfactant domains are subject to 
coarsening. The mobility of the surfactants during this 
process is dependent on the height of the liquid film un- 
derneath (see equation (J2|). We start with a systematic 
investigation of the synchronization behavior of the one 
dimensional system. To isolate the effect of the prepat- 
terned substrate on the time-periodic formation of sur- 
factant domains from the coarsening, we first determine 
the location Xf behind the contact line, where the do- 
mains are created (see figure [3|. Note, that this loca- 
tion depends on the transfer velocity V, since the menis- 
cus is more and more elongated with growing V. Then 
s(T) := T(Xf,T) is sampled as a time series. By use 
of a standard procedure from signal analysis, one can 
unambiguously define an instantaneous phase (j) s (T) of 
the "signal" s(T). To this end, one constructs the ana- 
lytic signal C(T) = s(T) +is H (T) =: C(T) exp(i0 a (T)), 
where sjj(T) is the Hilbert transform of s(T) [24]. Fi- 
nally, we obtain the frequency of the domain formation 
as lj = (d(/> s /dT), where (...) denotes the time average. 
This frequency can then be compared to the time periodic 
forcing frequency Q = 2ttV/L p due to the prepatterned 
substrate moving underneath the position Xf. 
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FIG. 2. Spacetime plots of solutions on homogeneous sub- 
strates (p = 0) at different transfer velocities, obtained from 
ID simulations: (a) V = 0.385 (b) V = 0.500 (c) V = 0.700. 
Only the monolayer density F is shown. Blue (dark) indicates 
the liquid expanded phase, while yellow (bright) corresponds 
to liquid condensed parts of the monolayer. 
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FIG. 3. Patterns of equal wavelength, Ao ~ 200, obtained 
from ID simulations at different velocities V on a homoge- 
neous substrate: (a) V = 0.4 and (b) V = 0.795. The red 
graph (solid line) corresponds to monolayer density T, while 
the blue graph (dashed line) represent the height profile H. 
The patterns are almost inverse to each other. The solid ver- 
tical lines indicate the location Xf where new domains are 
perpetually created. 



FIG. 4. Wavenumber ko of patterns transfered onto homoge- 
neous substrates against transfer velocity V, as obtained from 
ID computations. The curve is nonmonotonous, so that the 
same ko can be obtained for high and for low V. The origins 
of the Arnold tongues in figure [7] can be determined from the 
intersection points of the curve with horizontals at the corre- 
sponding n : m ratio. This is shown here for 1 : 1 and 3 : 2 
synchronization with L p — 200. 



B. Choice of parameters 



The substrate-liquid interaction is quantified by the 
Hamaker constant as. Here, we use a typical value 
as = 0.5/(127r) • 10 -20 J and assume a precursor height 
of hp « 5.8 nm, yielding the dimensionless Hamaker con- 
stant A = 6.5 • 10~ 2 . 

The critical number density of the first order phase 
transition between the LE and the LC phase is taken to 
be 7 cr = 1.539 • 10 18 ra -2 . We assume the surface ten- 
sion of water covered with a monolayer of density j CY to 
be 10% lower than the surface tension in absence of any 
surfactant: ao = 0.9a a b s = 65.48 mN/m. Then, the ex- 
perimental lateral pressure isotherm obtained for DPPC 
at 20 ° C is reasonably well approximated by using the 
parameters Mi = 50 and = 0.2 in the equations ^ 
and (16). Estimating the surfactant domain line tension 
k = 5 • 10 -52 Jm 4 , we obtain the time and length scales 
lo = 134.49 nm and to = 25.23 /is. This yields the small- 
ness parameter e = 4.31 • 10 -2 . 

The evaporation parameters used in our simulation are 
given by Ev = 0.45 and /i v = 5 • 10 -4 , amounting to a 
rate constant Q e = 4.93 • 10 -9 m 2 s/kg. 

The SMC coupling constant is set to B = 1500, corre- 
sponding to a lateral pressure shift Ap\ at = S(H p ) ao = 
-3ABa /10 « -3.60 mN/m caused by SMC. Thus, the 
coexistence lateral pressure at the substrate is roughly 
5% lower than for the floating monolayer. 



C. Homogeneous substrates 

On a homogeneous substrate, p = 0, spatiotemporal 
oscillations are found in a velocity range 0.38 < V < 0.83. 
Since the corresponding wavenumbers and frequencies 



are intrinsic quantities of the system, they are denoted 
by fco an d cjo, respectively, in order to distinguish them 
from the wavenumbers k and frequencies uj obtained on 
prepatterned substrates. Figure [4] shows ko = uio/V as 
a function of transfer velocity, obtained from one dimen- 
sional calculations. Note, that ko(V) is non- monotonous, 
yielding pairs of patterns with equal k obtained at differ- 
ent velocities. The two patterns comprising such a pair 
are roughly inverse to each other, but differ in detail. 
This is shown in figure [3] where two patterns of equal 
wavelength, Ao ~ 200, exhibit different structure: The 
LE domains in figure [3] (a) are narrower than the LC do- 
mains in figure [3] (b). As will be shown in the next para- 
graph, these different types of patterns will react slightly 
different to the periodic forcing caused by prepatterned 
substrates. 



For most velocities, the two dimensional system be- 
haves analogously, yielding stripes parallel to the contact 
line, corresponding to the alternating domains of the one- 
dimensional case. These solutions are symmetric with re- 
spect to translations along the F-axis, that is, the system 
behaves just as expected from the results obtained in the 
one dimensional case. However, for values of V chosen 
close to the upper or lower bound of the velocity range 
where patterns are obtained, the stripe formation is not 
stable. Instead, after a few stripes are produced, the 
regularity is destroyed and disordered domains are cre- 
ated. Remarkably, this disordered state is only transient 
for velocities chosen from the range V ^ 0.47 and marks 
an orientation transition from parallel to perpendicular 
stripes [T4] . 
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FIG. 5. Synchronization plateaus obtained from ID simula- 
tions of transfer onto a prepatterned substrate with L p — 400 
and p = 0.006. 
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FIG. 6. Staircase obtained from one dimensional simulations 
using L p = 400 and p = 0.006, after elimination of the transfer 
velocity in favor of the natural frequency scale Q/cjo- Red 
crosses indicate patterns obtained for velocities higher than 
0.52, black squares correspond to lower V . 



D. Prepatterned substrates 

For p > the wettability contrast of the moving sub- 
strate leads to an oscillation of the meniscus, going along 
with a change of the wavelength of the produced pattern. 
Within certain velocity ranges, the pattern synchronizes 
with the substrate, yielding commensurate frequency ra- 
tios uj/Q = n : m with integer n, m. This can be seen 
from the pronounced plateaus shown in figure [5j It is 
also noteworthy, that the overall velocity range where 
patterns are obtained is extended beyond V = 0.83, the 
upper limit on homogeneous substrates. As an alterna- 
tive to the representation used in figure |5j one can also 
use the ratio Q/uo as the abscissa, thereby directly re- 
lating the observed frequencies to the natural frequencies 
of the system without referring to the transfer velocity 
V. In the case considered here, this approach has two 
disadvantages. Firstly, velocities V > 0.83, where pat- 
terns are only obtained on prepatterned substrates, could 
not be included in the diagrams, since for them ujq = 0. 



Secondly, as has been mentioned in explanation of fig- 
ures [4] and |3J patterns obtained for higher V differ from 
those obtained for lower V, even if they have exactly the 
same period. Due to these differences, the synchroniza- 
tion plateaus of the low and the high transfer velocity 
regimes do not match perfectly. Instead, one obtains two 
staircases, one for low and one for high V, as is shown 
in figure [6] Therefore, we will keep the V representa- 
tion throughout the remainder of this article. Naturally, 
the frequency plateaus broaden with increasing ampli- 
tude of the wettability contrast p. This dependence can 
be visualized in a p-V diagram, where the synchroniza- 
tion regions are drawn for different values of p, yielding 
the famous Arnold tongues [24] . Figure [7] shows these 
diagrams for the three different considered wavelengths 
L p . Generally, an n : m Arnold tongue meets the V- 
axis at the velocity corresponding to the n : m ratio of 
the intrinsic wavenumber ko and the forcing wavenumber 
kQ. This velocity can be easily found by looking at the 
intersections of a horizontal line at ko = k^n/m with the 
curve ko(V). This is shown in figure [4] for the examples 
of 1 : 1 and 3 : 2 synchronization with L p = 200. In 
comparison to the case of homogeneous substrates (see 
figures [2] and [3| , transfer onto prepatterned substrates 
yields patterns of higher complexity. Figure [8] shows 
snapshots of transfer using a periodic wettability con- 
trast of wavelength L p = 400 and amplitude p = 0.005. 
One has to distinguish two cases. For velocities which 
belong to one of the synchronization plateaus, the ob- 
tained pattern is periodic. Figure [8] a) for example shows 
a periodic train of domain pairs produced by 2 : 1 syn- 
chronization. If the chosen velocity instead lies outside 
the plateaus, a nonperiodic train of domains is obtained 
(see figure [8] b)). By eye, it is generally difficult to dis- 
tinguish a non-synchronous pattern from a synchronous 
pattern with higher values of n and m as can be seen by 
comparison of figures [8] b) and[8]c). For the system un- 
der consideration, this difficulty is even enhanced since, 
as has been mentioned above, the domains are subject to 
coarsening from the moment they detach from the con- 
tact line. The coarsening is obvious in figure [8] d), where 
the leftmost, that is oldest, pair of domains already dif- 
fers visibly from the rightmost pair which has just been 
formed at the meniscus. 

Like on homogeneous substrates, in the case of two 
spatial dimensions, stripe patterns arranged parallel to 
the contact line are obtained analogously to the domains 
found in the one dimensional case. The periodicity of 
these patterns matches the one dimensional solutions. 
Figure [9] shows two snapshots of numerical simulations 
for a prepattern of period L p , However, as has been 
mentioned above in the discussion of transfer onto homo- 
geneous substrates, there exists a velocity range, where 
the formation of stripes arranged parallel to the contact 
line competes with a stationary pattern of perpendicular 
stripes. Due to the presence of the prepattern, this com- 
petition gives rise to the formation of new intermediate 
patterns, which are neither stripes parallel nor stripes 
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FIG. 7. Arnold tongues resulting from one dimensional com- 
putations for prepatterned substrate with (a) L p = 200, (b) 
L p = 300, and (c) L p = 400. The dashed lines show how the 
Arnold tongues are expected to go to p = 0. The vertical solid 
line indicates the velocity corresponding to the maximum of 
the intrinsic frequency ojq. 



1 



(a) 2:1, V=0.405 







500 



1000 
X 



1500 



2000 



(b) 



V=0.685 



1 JWzJLAi\^i^^ 



500 



1000 
X 



1500 



2000 



u 



(c) 5:2, V=0.72 



500 



1000 
X 



1500 



2000 



u 



(d) 2:1, V=0.815 



500 



1000 
X 



1500 



2000 




2000 



FIG. 8. Patterns in the ID system, obtained during transfer 
with five different velocities on a prepatterned substrate with 
L p = 400, p — 0.005. The velocities in (a) and (c-f) are chosen 
from different n : m synchronization regions, while (b) shows 
a pattern obtained at a velocity from between the right 3 : 1 
and 2 : 1 regions (compare [7] (a)). 
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FIG. 9. Two-dimensional patterns obtained on a prepatterned 
substrate with L p = 400 for two different transfer velocities 
V. The height information corresponds to the film profile H, 
while the monolayer density T is color-coded. 
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FIG. 10. Two dimensional solutions in the velocity range of 
the parallel-perpendicular orientation transition for various 
wetting contrasts p. 

perpendicular to the contact line and which are not ob- 
served on homogeneous substrates. This is explained in 
detail in the next subsection. 



E. The orientation transition 

During transfer onto homogeneous substrates, a transi- 
tion from stripes parallel to the contact line to perpendic- 
ular stripes was observed for velocities 0.38 < V < 0.47. 
On prepatterned substrates, these perpendicular stripes 
are stable only for small wetting contrast p. To illustrate 
this, we consider the solution obtained for p = 0, V = 
0.44, that is a stable, stationary pattern of stripes per- 
pendicular to the contact line obtained on a homoge- 
neous substrate. This solution is then used as the initial 
condition for computations with non-zero p. In these 
simulations, the wetting contrast is not activated instan- 
taneously across the whole domain. Instead, the region 
with p > is moving in from the right boundary with 
the transfer velocity V, thereby mimicking transfer onto 
a substrate, which is divided into an homogeneous and 
a prepatterned part. This means, that the monolayer is 
first transfered onto a homogeneous substrate, where it 
forms stripes perpendicular to the contact line. Then, at 



T = 0, the contact line reaches the prepatterned part of 
the substrate. Figure [To] shows, how the wetting contrast 
impacts the initially regular stripe pattern. For small 
amplitudes such as p = 0.002, the stripes persist but 
are periodically undulated across the X-direction. These 
undulations travel along with the moving prepattern un- 
derneath. For higher values of p the stripes break up 
into well- arranged arrays of circular LE domains, aligned 
like beads on a string parallel to the contact line, and 
advected with the substrate. For even higher amplitudes 
p ~ 0.007, these linearly arranged LE domains alternate 
with solid LE stripes. Finally, for very strong wetting 
contrasts p > 0.01, the extrinsic structure of the prepat- 
tern is completely enforced upon the system, resulting in 
a regular L p -periodic stripe pattern . 



IV. CONCLUSIONS AND OUTLOOK 

A theoretical investigation of monolayer transfer onto 
chemically prepatterned substrates has been presented. 
In our model, a spatially varying disjoining pressure ac- 
counts for the wetting heterogeneity due to the prepat- 
tern. Focussing on the case, where the wetting properties 
change periodically in the direction perpendicular to the 
contact line, complex periodic patterns were obtained as 
a result of synchronization effects of various orders. In 
addition, it has been shown, that a sufficiently strong 
wetting contrast can be used to inhibit the transition 
from stripes parallel to the contact line to perpendicu- 
lar stripes. Also, for moderate wetting contrasts, well- 
ordered arrays of circular LE domains are obtained, as 
well as mixed patterns consisting of alternating stripes 
and circular domains. These results demonstrate, that 
even simple prepatterns provide means to create vari- 
ous complex, well-ordered patterns in monolayer transfer 
systems and therefore open new perspectives for the con- 
trolled nanopatterning of surfaces. 

It can be expected, from our findings, that patterns of 
higher complexity might be obtained by use of more com- 
plicated prepatterns which are not necessarily structured 
periodically. But even using a periodic wetting hetero- 
geneity as is considered in this article, interesting new 
patterns might emerge, if the substrate is withdrawn un- 
der another angle, so that the stripes of the prepattern 
are no longer parallel to the contact line. These possibili- 
ties remain to be elaborated in a joint effort of theoretical 
and experimental research. 
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